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We convey an idea of the significant recent progress, which opens up good perspectives for 
high-precision ab-initio computations in heavy flavour physics based on lattice QCD. Rather 
than surveying the latest results, this contribution focuses on the concept and the challenges of 
fully non-perturbative computations in the B-meson sector, where the b-quark is treated within 
an effective theory. We outline its use to determine the b-quark mass and report on the results 
obtained in the quenched approximation and on the status in the two dynamical flavour theory. 

1 B-physics and lattice QCD 

The plenty of beautiful results from recent and still ongoing B-physics experiment^El 
which require the knowledge of QCD matrix elements for their interpretation in terms of 
parameters of the Standard Model and its possible extensions, motivates investigations in 
lattice QCD. The importance of this interplay of experiment and theory is further expressed 
by the fact that one of its main objectives, the phenomenon of CP violation, is closely 
related to the symmetry breaking mechanism the 2008's Nobel Prize was dedicated to. 

Lattice QCD represents our best founded theoretical formulation of QCD and allows 
for the computation of low-energy hadronic properties in the non-perturbative domain, 
where the usual power series expansion in the coupling constant fails, through the Monte 
Carlo evaluation of the Euclidean path integral after a discretization of space-time on a 
lattice with spacing a in all 3+1 dimension^Ml, While such numerical computations nec- 
essarily involve approximations, one of the key features of the lattice approach is that all 
approximations can be systematically improved. For an overview of results from the field 
of heavy flavour physics, which reflect some of these improvements by the small error bars 
quoted for many quantities, we refer to the reviews of past Lattice ConferenceJSEl. 

1.1 Challenges 

Among the various considerable challenges one faces in an actual lattice QCD calculation 
on the theoretical and technical levels, let us only highlight the multi-scale problem, which 
is also particularly relevant in view of B-physics applications. This is illustrated in Fig.[T] 
There are many disparate physical scales to be covered simultaneously, ranging from the 
lightest hadron mass of TOtt ~ 140 MeV over tod ~ 2 GeV to me ~ 5 GeV, plus the 
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Figure 1 . Large range of energy (p) scales in lattice QCD, where shaded areas refer to quark mass values (in the 
MS scheme) quoted by the Particle Data Group**. Red marks indicate the pion, the D- and the B-meson mass. 

ultraviolet cutoff of Auv = a^^ of the lattice discretization that has to be large compared 
to all physical energy scales for the discretized theory to be an approximation to the con- 
tinuum one. Moreover, the finiteness of the linear extent of space-time, L, in a numerical 
treatment entails an infrared cutoff Air = L^^ so that the following scale hierarchy is met: 

Air = L^^ <c m^r , . ■ . , tod , tob < a^^ = Auv • (1) 

This implies L > A/m-,, ss 6 fm to suppress finite-size effects in the light quark sector and 
a < 1/ (2mD) ~ 0.05 fm to still properly resolve the propagation of a c-quark in the heavy 
sector Lattices with L/a > 120 sites in each direction would thus be needed to satisfy 
these constraints, and since the scale of hadrons with b-quarks was not even included to 
arrive at this figure, it is obvious that the b-quark mass scale has to be separated from the 
others in a theoretically sound way before simulating the theory. In Sec. |2] we describe, 
how this is achieved by recoursing to an effective theory for the b-quark. 

Another non-trivial task is the renormalization of QCD operators composed of quark 
and gluon fields, which appear in the effective weak Hamiltonian, valid at energies far 
below the electroweak scale. Besides perturbation theory^, powerful non-perturbative ap- 
proaches have been developecffil, and we will come back to the non-perturbative subtrac- 
tion of power-law divergences in the context of the effective theory for the b-quark later. 

1.2 Perspectives 

As for the challenges with light quarks, we only mention that the condition L > 6 fm may 
be relaxed by simulating at unphysically large pion masses, combined with a subsequent 
extrapolation guided by chiral perturbation theory'" and its lattice-specific refinements. 

Regarding the algorithmic side of a lattice QCD simulation, the Hybrid Monte CarlcP 
(HMC) as the first exact and still state-of-the-art algo rithm has received considerable im- 
provements by multiple time-scale integration scheme^-^Ii^ the Hasenbusch trick of mass- 
preconditionin^nnH, supplemented by a sensible tuning of the algorithm's parameter^, 
and the method of domain decomposition (DD) applied to QCD'^^IISI^ just to name a few. 
In addition, low-mode deflation^" (together with chronological inverteri^ has led to a 
substantial reduction of the critical slowing down with the quark mass in the DD-HMC. 

Finally, in parallel to the continuous increase of computer speed (at an exponential 
rate) over the last 25 years and the recent investments into high performance computing 
at many places of the world, the Coordinated Lattice Simulations^^ (CLS) initiative is a 
community effort to bring together the human and computer resources of several teams 
in Europe interested in lattice QCD. The present goal are large-volume simulations with 
iVf = 2 dynamical quarks, using the rather simple 0(a) improved Wilson action"* to profit 
from the above algorithmic developments such as DD-HMC, and lattice spacings a = 
0.08, 0.06, 0.04 fm, sizes i = (2 - 4) fm and pion masses down to = 200 MeV, 
which altogether help to diminish systematic and statistical errors. Amongst others, charm 
physic^ as well as our B-physics programme outlined here are being investigated. 
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2 Non-perturbative Heavy Quark Effective Theory 



Heavy Quark Effective Theory (HQET) at zero velocity on the lattice offers a reliable 
solution to the problem of dealing with the two disparate intrinsic scales encountered in 
heavy-light systems involving the b-quark, i.e., the lattice spacing a, which has to be much 
smaller than 1 / rrib to allow for a fine enough resolution of the states in question, and the 
linear extent L of the lattice volume, which has to be large enough for finite-size effects to 
be under control (recall also Fig.[T]l. 

Since the heavy quark mass (rrib) is much larger than the other scales such as its 3- 
momentum or Aqcd ~ 500 MeV, HQET relies upon a systematic expansion of the QCD 
action and correlation functions in inverse powers of the heavy quark mass around the 
static limit (mb ~^ oo). The lattice HQET action at 0(l/mb) reads: 

S'hqet = a^Y.x'^i-i {Do + 5m- WkinD^ - cjspinCrB} V'h , (2) 

with -0h satisfying P+V'h ~ V'h, P+ = ^^2° ' parameters Wkin and Wspin being 

formally 0(l/mb). At leading order (static limit), where the heavy quark acts only as 
a static colour source and the light quarks are independent of the heavy quark's flavour 
and spin, the theory is expected to have ^ 10% precision, while this reduces to ~ 1% at 
0(l/TOb) representing the interactions due to the motion and the spin of the heavy quark. 
As crucial advantage (e.g., over NRQCD), HQET treats the 1/mb-corrections to the static 
theory as space-time insertions in correlations functions. For correlation functions of some 
multi-local fields O and up to l/rrib-corrections to the operator itself (irrelevant when 
spectral quantities are considered), this means 

(O) = (O)stat + 51 Kin(OOkin(x))stat + C^spin (©apin (2:)).tat } , 0) 

X 

where (C')stat denotes the expectation value in the static approximation and Okin and 
Cspin are given by V'hD^V'h and ilj^a'BipY^. In this way, HQET at a given order is (power- 
counting) renormalizable and its continuum limit well defined, once the mass counterterm 
6m and the coefficients Wkin and Wspin are fixed non-perturbatively by a matching to QCD. 

Still, for lattice HQET and its numerical applications to lead to precise results with 
controlled systematic errors in practice, two shortcomings had to be left behind first. 

1. ) The exponential growth of the noise-to-signal ratio in static-light correlators, which 
is overcome by a clever modification of the Eichten-Hill discretization of the static actiorP^. 

2. ) As in HQET mixings among operators of different dimensions occur, the power- 
divergent additive mass renormalization 5m ~ On/a already affects its leading order. 
Unless HQET is renormalized non-perturbatively^, this divergence — and further ones 
^ g^jc? arising at 0(l/mb) — imply that the continuum limit does not exist owing to 
a remainder, which, at any finite perturbative ordei^^ diverges as a — >^ 0. A general 
solution to this theoretically serious problem was worked out and implemented for a deter- 
mination of the b-quark's mass in the static and quenched approximations as a test casd^ll 
It is based on a non-perturbative matching of HQET and QCD infinite volume. 

3 Application: The b-quark mass from HQET at 0(1 /mb) 

Let us first notd^ that in order not to spoil the asymptotic convergence of the series, the 
matching must be done non-perturbatively — at least for the leading, static piece — as 
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soon as the 1 /mb-corrections are included, since as mb oo the perturbative truncation 
error from the matching coefficient of the static term becomes much larger than the power 
corrections ^ Aqcd /^rib of the HQET expansion. 

QCD HQET 




Figure 2. Idea of lattice HQET computations via a non-perturbative determination of HQET parameters from 
small- volume QCD simulations. For each fixed Li, the steps are repeated at smaller a to reach a continuum limit. 



In the framework introduced in , matching and renormalization are performed simul- 
taneously and non-perturbatively. The general strategy, illustrated in Fig. |2j can be ex- 
plained as follows. Starting from a finite volume with Li w 0.5 fm, one chooses lattice 
spacings a sufficiently smaller than 1 /mb such that the b-quark propagates correctly up 
to controllable discretization errors of order a^. The relation between the renormalization 
group invariant (RGI) and the bare mass in QCD being known, suitable finite-volume ob- 
servables $^(^1, A/h) can be calculated as a function of the RGI heavy quark mass, M\y, 
and extrapolated to the continuum limit. Next, the power-divergent subtractions are per- 
formed non-perturbatively by a set of matching conditions, in which the results obtained 
for are equated to their representation in HQET (rh.s. of Fig.[2]i. At the same physical 
value of Li but for resolutions Li/a = 0(10), the previously computed heavy-quark mass 
dependence of $fe(Li, M^) in finite-volume QCD may be exploited to determine the bare 
parameters of HQET for a sa (0.025 — 0.05) fm. To evolve the HQET observables to large 
volumes, where contact with some physical input from experiment can be made, one also 
computes them at these lattice spacings in a larger volume, L2 = 2Li. The resulting re- 
lation between and ^k{L2) is encoded in associated step scaling functions (SSFs) 
(Tfc, indicated in Fig.|2] By using the knowledge of ^kiL2, M^) one fixes the bare parame- 
ters of the effective theory for a « (0.05 — 0.1) fm so that a connection to lattice spacings is 
established, where large-volume observables, such as the B-meson mass or decay constant. 
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can be calculated (bottom of Fig. [2]). This sequence of steps yields an expression of niB, 
the physical input, as a function of A/h via the quark mass dependence of <i>/j(ii, M^), 
which eventually is inverted to arrive at the desired value of the RGI b-mass within HQET. 
The whole construction is such that the continuum limit can be taken for all pieces. 



3.1 Review of the quenched computation of the b-quark masd^ 

To apply this to A/b, the task is to fix Sm and oj^in non-perturbatively by performing a 
matching to QCD, after restricting to spin-averaged quantities to get rid of the contributions 
proportional to Wspin, For sensible definitions of the required matching observables, $i and 
$2, we work with the Schrodinger functional (SF), i.e., QCD with Dirichlet boundary con- 
ditions in time and periodic ones in space (up to a phase for the fermions): ^^^'^ [L, mh) 
exploits the sensitivity of SF correlation functions to 6 and $^'~^°(L, mh) = Lri(L, mh), 
where Fi is a finite- volume effective energy. When expanded in HQEfl $f*^'^'^(L) 
is given by Wkin times a quantity defined in the effective theory (called 6*, 6*')), 

whereas $^^^^(L) is a function of Wkin and mbaic — Sm + rrih involving two other 
HQET quantities, rf'**(L) and T\"^{L). According to the strategy sketched above, by 
equating $|?^^(Li,mh) and $^^^"^(Li) one can determine the bare parameters mbaro 
and Wkin as functions of mh at the lattice spacings belonging to the volume Lf. To employ 
the spin-averaged B-meson mass, mg', as phenomenological input, the are evolved to 
larger volumes through proper SSFs, where the resulting $^'^^^(2Li, mh) still carry the 
dependence on mh inherited from the matching to QCD in Lf. After 2 evolution steps 
(and taking continuum limits), linear extents of > 1.5 fm are reached, and mbarc and Wkin, 
expressed in terms of SSFs, $^^^(Li, mh) as well as i?^'"' ^f^^ ^^'^ r!^"\ are obtained 
— again as functions of mh. Now, the b-quark mass is extracted by solving 

TO- ^ + Wkin(mh)S'^'" + mbare(mh) (4) 

for mh, with S^t^t = limi^oo Tf ^* and E^"' - - ( B | Oki„(0, z) \ B )stat. All 

quantities entering eq. Q have a continuum limit either in QCD or HQET, which implies 
that all power divergences have been subtracted non-perturbatively. 

In case of the leading-order, static approximation, where only mbarc needs to be de- 
termined, the small- and large-volume matching conditions simplify to ri(ii,mh) = 
rf^^{Li) + mbarc and m|^ — E^^^^ + mbarc, respectively. To be able to solve the first 
equation for mbare and replace it in the second, we bridge the volume gap in two steps by 
inserting a SSF a^lLi) = 2Li [Ff '^*(2Li) — Ff ^*(Li)] and arrive at the master equation 

Li [mg^ - {E-'-' - Ff-*) ] - = LiFi(Li, mh) , (5) 

where Fi originates from QCD in i| and any reference to bare parameters has disappeared. 
Its graphical solution is reproduced in Fig.|3]and yields A/^'*^' = 6.806(79) GeV. 

The inclusion of the sub-leading 1/mb-effects is technically more involved and ex- 
ploits the freedom of choices for the angle(s) 6 and an alternative set of matching observ- 
ablePl. We just quote the final value m|^^(mb) = 4.347(48) GeV with the remark that, 
upon including the 1 /mb-terms, differences among the static results w.r.t. the matching 
condition chosen are gone, which signals practically negligible higher-order corrections. 



"Here, Sm = in the action; its effect is accounted for in tlie overall energy sliift rribarc in HQET versus QCD. 
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Figure 3. Graphical solution of eq. js} in tlie quenched casJSI The quantity used in the finite-volume matching 
step is (Li , M) = LiTi (Li , M), where z = LiM and M = Af^ the RGI heavy quark mass. 

3.2 Status in two-flavour QCD 

The renormalization of HQET through the non-perturbative matching to A^f = 2 QCD in 
finite volume, to do the power-divergent subtractions, is under waj^" SI As an important 
prerequisite, the calculated non-perturbative relation between the RGI and subtracted bare 
heavy quark mas J23! enables to fix RGI heavy quark masses in the matching volume Lf: 

LiM = 2'M(.go)^(.go) (1 + V(.go)an^q,h) ii"^q,h , M = Mh . (6) 

The extent Li is defined via a constant SF coupling, g'^{Li/2) ~ 2.989, and the PCAC 
masses of the dynamical light quarks are tuned to zero. 

Fig.|4]shows two examples for the heavy quark mass dependence of finite-volume QCD 
observables in the continuum limit, which enter the non-perturbative matching. 
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Figure 4. ^-dependence of two QCD observables, which enter the matching of QCD with HQET in A^f = 2 with 
massless dynamical quarks, z = LiAf as before. Left: Spin-averaged B-meson energy. Right: Ratio of axial- 
vector to vector current matrix elements. The conversion functions Cx, which only introduce small perturbative 
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translate the HQET predictions to the corresponding QCD quantities at finite values of z. 



These results also allow to perform non-perturbative tests of HQET in the spirit of the 
corresponding quenched investigation^^. The calculation of the step scaling functions in 
HQET is expected to be finished soon, and the large-volume part of our strategy is currently 
being implemented within the CLS effori^Sl. 
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4 Outlook 



The non-perturbative treatment of HQET including l/wb-terms can lead to results with 
unprecedented precision for B-physics on the lattice. It also greatly improves our confi- 
dence in the use of the effective theory. The striking agreement, for example, between the 
decay constant Fb, computed including 1 /mb-corrections and the value resulting from 
the interpolation between the static number and data around the charnl^, though still in 
the quenched approximation, provides a strong internal check of the approach. In addition, 
the HQET parameters at 0(l/mb) calculated non-perturbatively by the ALPHA Collabo- 
ratiorJSlcan be employed for several other quantities. The programme aiming to reach the 
same accuracy in the iVf = 2 dynamical case is well advanced and progressing fast^*^. 
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